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The mechanism underlying the substantial amplification of the high-order harmonics q ± 2K (K 
integer) upon the addition of a weak seed XUV field of harmonic frequency qui to a strong IR field 
of frequency ui is analyzed in the framework of the quantum-mechanical Floquet formalism and 
the semiclassical re-collision model. According to the Floquet analysis, the high-frequency field 
induces transitions between several Floquet states and leads to the appearance of new dipole cross 
terms. The semiclassical re-collision model suggests that the origin of the enhancement lies in the 
time-dependent modulation of the ground electronic state induced by the XUV field. 

PACS numbers: 



03.65.-w, 42.50.Hz, 42.65.-Ky, 32.80.Rm 

Focusing intense linearly-polarized monochromatic IR laser pulses into gas of atoms can lead to the emission of high- 
energy photons with frequencies extending into the extreme ultraviolet (XUV) and X-ray region by high harmonic 
generation (HHG). The HHG phenomena stands as one of the most promising methods of producing short attosecond 
pulses (as-pulses) 1|. 



The contamination of the strong IR field with a second 



weak XUV fields has a dramatic 



effect on the dynamical behavior of the electrons, and had drawn a lot of attention in recent years. On the basis of the 



three-step (re-collision) model 8|, |9|, [l0| it had been argued that the role of the XUV field is to switch the initial step 
in the generation of high harmonics from tunnel ionization to the more efficient single XUV-photon ionization. This 
might explain the improved macroscopic HHG signal obtained in experiments: the XUV-assisted ionization increases 
the number of atoms which participate in the HHG process and improves phase matching 11 1 . 



The effect at the single-atom level, however, is less clear. It has been shown that the XUV photons control the 
timing of ionization, and preferentially select certain quantum paths of the electron 12| . While this effect may lead to 



the enhancement of the low-order harmonics in the plateau, it can't account for the large enhancement in the cutoff 



and beyond (FigQJ. A 3-step-model classical analysis of HHG suggests that the contribution of the XUV field to the 
kinetic energy of the returning electron is negligible. The kinetic energy of a classical free electron of charge e and 
mass m, driven by a linearly-polarized strong IR fundamental field of frequency u>, amplitude an d polarization 
e k (Ei(f) = eke 1 / 1 cos(u)i)) is Ek(t) = ^J^-. pit) = ^-[sin{ujt) — sin(u>ti)] is the momentum of the electron, and 
it has been assumed that the electron is freed at time Aj with zero momentum. The addition of a weak harmonic 
XUV field of frequency qui (where q is a large integer) and amplitude (e™ << e™) with the same polarization, 
(E q (f) = ek£g" cos(qujt)) , adds a small correction to the momentum, which is proportional to As a result, the 
correction to the kinetic energy, which appears in the form of two additional terms, proportional to and (-^j) 2 , 
is negligible. Thus, the additional XUV field will not affect the electron trajectories and will not contribute to their 
kinetic energy. For this reason the relative phase between the two fields doesn't play a role in the HGS, which is indeed 
verified in both classical analysis and quantum mechanical simulations (a small q, however, will affect the dynamics 



differently 2J, |f3|). In addition, assigning the electron a non-zero initial momentum to account for the photoelectric 
effect, will not increase its kinetic energy upon recombination. 

An illustrative TDSE simulation however (FigfTJ) shows an enhancement of the cut-off harmonics and the harmonics 
q ± 2K (K integer) upon addition of a weak XUV field to the strong IR field. Moreover, the HGS possesses certain 
symmetries: with respect to its center at harmonic q, the distribution of harmonics of the enhanced part of the 
spectrum (harmonics that have been produced only due to the addition of the XUV field) , is symmetric with respect 
to q and remains almost invariant upon variation of q. This suggests that despite the fact that the additional weak 
XUV field doesn't affect the electron trajectories, it does affect the recombination process. As will be shown later, 
the XUV field induces periodic modulations to the remaining ground electronic state, with the same frequency as the 
XUV field. The returning electronic wavepacket recombines with this modulated ground state to emit new harmonics. 
The purpose of this article is to reveal this mechanism which is responsible for the amplification phenomena due to 
the inclusion of the weak XUV field and to prove that the enhancement is a robust single atom phenomenon. The 
mechanism could suggest new types of HHG experiments. It is not limited to the description of the self-occurring 
case in monochromatic HHG experiments where XUV radiation, generated by the leading edge of the IR pulse, co- 
propagates with the IR field to form a bichromatic driver field in the last part of the medium (thus leading to the 
extension of the cutoff energy in real experiments as compared to single- atom calculations). In cases where the XUV 
field saturates, it might be useful to add it externally For example, He is known to produce higher harmonics than 



Ar. Hence, the support of the HGS obtained from Ar can be dramatically extended by shining the Ar with a high 
harmonic obtained from He (which is absent in the Ar HGS) in addition to the strong IR field. 




FIG. 1: (color online) HGS obtained from a ID model Hamiltonian of Xe atom irradiated by a 50-oscillation sine-square pulse of 
bichromatic laser field composed of a strong laser field of frequency uj (A = 800nm) and amplitude ej™ {If 1 — 4.3 • 10 13 W/cm 2 ) 
and a weak field of frequency qui and amplitude e," (7™ ~ 3.5 • 10 8 W/cm 2 ) for different values of q: q — 25 (solid red line), 
q — 37 (dotted green line), q = 52 (solid blue line). HGS in the absence of the XUV field is shown in the dotted black line 
where the position of the cutoff is at the 15th harmonic. The harmonics above the 29th harmonic, are enhanced in the addition 
of the XUV field, despite of its small intensity. In addition, with respect to its center q, the distribution of the new harmonics 
in the HGS is symmetric (i.e., for q — 37, a(33uj) ~ a(ilu>), etc.) and upon variation of q it shifts but remains almost invariant. 



In order to reveal the enhancement mechanism due to the inclusion of the weak XUV field, we study the dynamics of 
a single active electron in an atom described by the field-free Hamiltonian Ho(r) subjected to a long pulse of the IR field 
Ei(t) in the length gauge and under the dipole approximation. The long pulse evolves the system adiabatically 14 1 
from the initial ground state of the field-free Hamiltonian |</>i(r)) to a single resonance Floquet eigenstate (r,t)) 



of EqU] 



15j which describes the entire dynamics of the system. A formalism of time independent perturbation theory 



is applicable since the time t may be treated as an additional coordinate 
to, to', n, n' , M, K denote integers (g Z). 



1G 



1 71 ] . In the following, all the parameters 



Hf{v,t)\^ n {v,t)) = e$l\^%(r,t)) (1) 
where Hp(r, t) = Ho(r) — er ■ Ei(t) — ih^ is the Floquet Hamiltonian. The indices [j, to) label the eigenstates j 



within any given Brillouin zone m and r describes the internal degrees of freedom. The Floquct eigenfunctions of this 

19 1 (written in the usual dirac notation) (ipv^r, t)\ip^ m ,(r, t)) r ,t = 



operator satisfy the c-product inner product 



18 



5jj>5 mm i and form a complete set. Floquet eigenfunctions which lie within the m-th Brillouin zone may be defined 
as |V'j>„M)) = \i>f\r,t))e lulmt and (^M)| = (ipf ] (r,t)\e- iumt with energies sf ] m = sf ] +mhuj. The ket and 
bra Floquet eigenfunctions are periodic with period T = 2tt/lo and can therefore be decomposed as a Fourier sum 
\tpf\r,t)) = J2„ l^°n( r )) e<u,nt and(^ 0) (r,t)| = E„(<^> ' {r)\e~ lulni '■ Note that the Fourier components of the bra 



state are not complex-conjugated 



18] 



In order to calculate the HGS one may assume the Larmor approximation [20( and analyze the time-dependent 
acceleration expectation value a^(t) = JU(V^(r, t)|r|^>i (r, i)) r which is proportional to the emitted field. The ac- 



celeration in energy space is given by the Fourier transform (Q) = A dtaf^ (t)e Exploring only frequencies 



which are integer multiples of u> [SI = Mlo ( M g Z)], and using the property ^ Jq dt e lunt = <5„.o, the expression 



obtained is a[°'(Mw) = — u> 2 AI 2 En(v4°n*( r )i r l Vi°n+M( r )} r - It can t> e shown to be non-vanishing only for integer odd 



values of M , which is a well known feature of monochromatic HHG 



21] 



Suppose the weak XUV field E q (t) is added. A new Floquet problem is obtained, which could be described by 
the Floquet Hamiltonian H^ EW (v,t) = H [ p\r,t) + V(r,t), where the additional term V(r,t) = -er • E q (i) could 
be treated as a perturbation. Time-independent lst-order perturbation theory may be used to get an approximate 
solution for the Floquet Hamiltonian Hp En (r, t) as 



\i>? EW (r,t)) = \^\r,t)) + Yl 4'' W (s)l#M))e Wt (2) 

(j',m')#(l,0) 

where the coefficients c[ ,m (q) are given by 



1 [q) ~ 2 e£ " ek '^ e w _ £ (o) : mfRu ■ P) 

Using this solution the time dependent acceleration expectation value a.^ EW (t) = ^z(4>i EW ( r ) t)\ r \' l Pi EW (' r , t))r can 
be calculated. Keeping terms up to first order in e™ the following expression for the acceleration in the frequency 
domain is obtained: 



a™(M W )=a< 0) (M W ) 



(j',m')#(l,0) n 

This is the expression for the emitted HHG field. The HGS (a(Mui) = |af ew (Muj)\ 2 ) has the same features as 



those presented in FigfT] (see |22j|). The weak perturbative XUV field shifts the HGS beyond the cutoff obtained by 
the 1R field alone. In the Floquet formalism presented here the origin of the HHG enhancement phenomena lies in 
the interferences between the ground and excited Floquet states. The HGS is modified due to the dipole cross-terms 
introduced by the weak XUV field. 

The features in the HGS could also be explained in terms of the re-collision model. It was shown that the additional 
weak XUV field doesn't affect the electron trajectories, i.e., doesn't modify the kinetic energy of the re-colliding 
electron. According to the findings of the numerical simulation it must however affect the recombination process. 
To see this we turn into the semiclassical re-collision model 8| where the electronic wavefunction at the event of 
recombination could be described as a sum of the following continuum and bound parts. Under the strong field 
approximation the returning continuum part in the direction of the polarization ek (which we take as the x-direction 
from now on for simplicity) is a superposition of plane waves ij)\{x, t) = — 1= dktp c (k,t)e l ^ kx ^~i^ t ^ where k = fce x 
(k = |k|) is the momentum of the electron, Ek = 4^7" ^ s ^ ne USU& 1 dispersion relation and tp c (k,t) are expansion 
coefficients which weakly depend on time. It is assumed that the continuum wavepacket tp c (r,t) is separable in the 
x-coordinate and the 2 other lateral coordinates such that ^c(r, t) = ip\ (x, t)ip^r(y, z, t). It is assumed that the ground 
state is only slightly depleted during the tunnel-ionization and that due to the ac-Stark effect the electron adiabatically 
follows the instantaneous ground state of the potential which is periodically modified by the IR and XUV fields. 
Since the ac-Stark corrections to the instantaneous energy and wavefunction are small for normal field intensities, 
the instantaneous ground state could be approximated as ipb(r,t) = <j>\{x + e1 ut cos(ujt) + e° ut cos(qujt) 7 y, z)e + ^ Ipt 
(where I p > and <f>i(r) are the field-free ground state eigenvalue and eigenstate, respectively). Note that ^b(r, i) 
approximately describes the resonance Floquet state \i/j^ EW (r, t)). The quiver amplitudes e° ut ,e° ut of the spatial 
oscillations of the ground state are of the order of e° ut — -^j2, i.e., a tiny fraction of a Bohr radius for normal 
laser intensities and/or large values of q. The bound part may therefore be expanded in a Taylor serie as ipb{ r , t) — 
e+£ 7 p*{0 1 (r) + [e° ut cos(u)t) + e° ut cos(qu;t)]J^;(j)i(r)}. Using the total wavefunction at the event of recombination 
\l/(r, t) = ipb(r,t) + ip c (r,t), the time-dependent acceleration expectation value a(t) = ^(^(r, t)\ — VVb(r)|\f'(r, t)) r 
could be calculated, where Vb(r) is the field-free potential. The dominant terms that are responsible for the emission 
of radiation at frequencies other than the incident frequencies u> and quj are the bound-continuum terms a(i) ~ 
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23?(-06(r, t)\ — VVb(r)|-0 c (r, t)) r . After some algebra it could be realized that the acceleration is composed of oscillating 
terms of the form 



a(*) ~ -23?^ / dk[i> IR {k)e-* [Ek+I ? ]t + ^ xuv {k)e° 1 ut e-* [Ek+I ? ±hu}]t + ^ xuv {k)e° q ut e-^ Ek+I " ±qhuj]t } (5) 

V27T J-oo 

where i> m {k) = d 3 r0i(r)VVb(r)^ (y, z, t)e lkx and 4>xuv( k ) = 

■^ipc(k) d 3 r d ^g^ VVq {t)i^>c (y, z, t)e ikx and the ± sign in each of the last two terms stands for summa- 
tion over two-terms each. The emitted field in a single re-collision event is composed of a continuum of these 
frequencies. 

It is therefore seen that despite of their small magnitude, the periodic time-dependent modulations to the ground 
electronic state induced by the XUV weak field of frequency quj are responsible for the appearance of the new harmonics 
around q in the HGS via recombination with the returning electronic wavepacket. Each electron trajectory (plane 
wave) with kinetic energy E^, recombines with the nucleus to emit, with equal probabilities, one of three possible 
photons with energies: I p + Ek, qhw + I p + Ek or qfiuj — (I p + Ek) ■ The HGS in the presence of the IR field alone 
hfl = I p + Ek is now shifted by the energy of the XUV photon Kqu>, and new harmonics are also formed, such that their 
distribution about the center q is symmetric. Also, with respect to the center q, the distribution of the XUV-formed 
harmonics, is invariant to a change in the energy of the XUV photon hquj, since these harmonics are "born" from the 
same set of electron trajectories which are characteristic of the IR field alone. When each single re-collision event is 
repeated every half cycle of the IR field, integer harmonics q ± 2K are obtained in the HGS. To see this, note that in 
two consecutive re-collision events at times t r and t r +T/2 the following symmetry holds: ip c {k 1 t r + T/2) — tp c {—k, t r ). 
Consequently, since Vq(t) and </>i(r) are symmetric functions for atoms (and a ^ r ^ is antisymmetric), the following 
symmetry holds ip IR (k, t r +T/2) = —ip IR (—k, t r ). The acceleration which results from the IR field therefore switches 
signs between subsequent re-collision events, which is the origin of the odd-selection rules. However, the behavior of 
the coefficients resulting from the addition of the XUV filed is different ip xuv (k,t r + T/2) = +ipxuv(~k,t r ). The 
acceleration which results from the additional XUV field doesn't switch signs between subsequent re-collision events 
and therefore yields even harmonics around q. 

The above suggestion could be verified by plotting the time- frequency distribution of high harmonics (Figj2|) obtained 
from the time-dependent acceleration expectation value whose spectra is given in Figfl] for q = 52. In accordance 
with the classical re-collision model, different harmonics are emitted repeatedly every half cycle, with the IR cut- 



off harmonic (the 15th harmonic) emitted at times ~ 0.2T + K ■ 0.5T. At those instants, also the 38th and 66th 
harmonics, which are produced by the most energetic IR trajectory, are emitted. Each electron trajectory in general, 
which under the IR field alone produces an harmonic O, generates upon the addition of the XUV field, two duplicated 
new harmonics with energies qhio + KQ, qfnu — Ml, and similar properties. For example the harmonics of orders 38 — 48 
and 56 — 66 have a "plateau" character (constant intensity), like the plateau harmonics 5 — 15. Moreover, as Eq[5] 
predicts, the intensity ratio of the enhanccd-plateau harmonics and the IR-plateau harmonics should be (for any of 
the values of q given in FigQ} (e^"') 2 = {-^yi) 2 — 1CU 10 , in agreement with the results of Fig[TJ 
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FIG. 2: (color online) Top view (red color- high intensity, blue color- low intensity) of the absolute square of the Gabor- 

„50t (t-to) 2 _snt 

tansformed acceleration expectation value ( J q a(t)e -r 2 e 1 , r = 0.1T) of the quantum mechanical simulation 
described in Fig[T]for q = 52, as function of to and fl. 



In conclusion, we have shown that the addition of a weak XUV harmonic field to a strong IR field leads to the 
extension of the cut-off in the HGS. The results of the quantum analytical expressions, quantum numerical simulations 
and classical arguments suggest that the enhancement is a single-atom phenomena. The seed XUV field modulates 
the ground state and affects the recombination process of all returning trajectories, and leads to the generation of new 
harmonics with structure well related to the HGS in the presence of the IR field alone. This amplification mechanism 
for the generation of high-order harmonics might be used to enhance the yield of harmonics in HHG experiments. 
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